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Abstract. We present the observation of the Tycho supernova remnant obtained with the EPIC and RGS in- 
struments onboard the XMM-Newton satellite. We compare images and azimuthally averaged radial profiles in 
emission lines from different elements (silicon and iron) and different transition lines of iron (Fe L and Fe K). 
While the Fexvn L line and Sixni K line images are globally spatially coincident, the Fe K emission clearly peaks 
at a smaller radius, indicating a higher temperature toward the reverse shock. This is qualitatively the profile 
expected when the reverse shock, after travelling through the outer power-law density profile, has entered the 
central plateau of the ejecta. The high energy continuum map has an overall smooth distribution, with a similar 
extent to the radio emission. Its radial profile peaks further out than the lines emission. Brighter and harder 
continuum regions are observed with a rough bipolar symmetry in the eastern and western edges. The spectral 
analysis of the southeastern knots supports spatial variations of the relative abundance of silicon and iron, which 
implies an incomplete mixing of the silicon and iron layers. 
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1. Introduction 

The Tycho supernova remnant (SNR) is considered as the 
prototype for Type la remnants (Baade 1945). The su- 
pernova was recorded in 1572 by Tycho Brahe, hence, the 
remnant is relatively young and still retains much infor- 
mation about the progenitor. Its X-ray emission is dom- 
inated by shocked ejecta. With an angular diameter of 8 
arcmin, Tycho SNR is well suited for investigating the spa- 
tial distribution of heavy elements in the ejecta. Hwang & 
Gotthelf (1997) presented deconvolved line images of the 
object using ASCA data. They found that the Fe K radial 
profile peaks at smaller radius than the other emission 
lines. This was interpreted as a possible iron stratification 
in the interior of the ejecta. 

With XMM-Newton observatory, it is now possible 
to perform true spatially resolved spectroscopy of Tycho 
SNR, allowing the mapping of the emission for individ- 
ual elements (and individual ions) and for the continuum 
at high energy. In this paper we present a first study of 
Tycho SNR with XMM-Newton. 
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2. Observation 

Tycho SNR was observed by XMM-Newton (Jansen et al. 
2001) on 2000 June 29 with the two X-ray instruments 
EPIC (European Photon Imaging Camera; Struder et al. 
2001; Turner et al. 2001) and RGS (Reflection Grating 
Spectrometer; den Herder et al. 2001). After selecting 
good time intervals (when the light curve in background 
regions is stable, i.e. uncontaminated by flaring events), 
the available exposure time was 12 ks for each EPIC 
MOS camera and 15 ks for the EPIC PN camera. Due 
to the high signal-to- noise (S/N) ratio of the RGS spec- 
tra of Tycho and the smaller sensitivity of the RGS in- 
struments to background variations, no time-selection was 
applied to the RGS data, resulting in a total of 45 ks of 
exposure time. The X-ray image of Tycho is entirely con- 
tained within the central CCD of each of the MOS cameras 
whereas in the PN camera it is spread over 6 CCDs, re- 
sulting in data gaps between each CCD. For simplicity, 
therefore, we have only used the EPIC MOS data in our 
spatial analysis. 

The three EPIC spectra (MOS 1, MOS 2 and PN), 
integrated over the entire remnant, are shown in Fig. 1, 
with the RGS spectrum obtained on a region smaller than 
Tycho (the FOV is shown in Fig. 2). The most prominent 
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lines in the EPIC spectra are the He a lines of magnesium, 
silicon, sulphur, argon, calcium, iron and an Fexvn line in 
the Fe L complex. With the RGS, the remnant has been 
detected up to the 5th order, however with a low S/N ratio 
from the third order. Emission lines from Sixni He a, 
Mgxi He a, Nex Ly a, Neix He a, two Fexvn lines (at 
15 A and 17 A) and O Viii Ly a have been identified. The 
structure seen both in EPIC and RGS spectra between the 
magnesium and neon lines (10 A and 12 A) is attributed 
to Fe L shell lines. Both RGS data sets were analyzed and 
are fully consistent with each other. 

3. Elemental structure in the ejecta 

The bright emission lines of heavy elements (silicon, sul- 
phur, iron) in the spectrum of the entire remnant originate 
from shocked material within the supernova ejecta. 

We have created images selecting photons in narrow 
energy bands corresponding to the main line features. For 
the maps of Sixin and Fexvn, we did not substract the 
continuum: these two lines have the highest S/N ratio 
above the continuum. For the lower S/N ratio Fe Ka line, 
where the spectral shape is steep, the continuum does need 
to be subtracted. This was done by using the overall global 
spectrum to derive the ratio of the continuum at the line 
energy to the continuum in an adjacent energy band, free 
of contamination by line emission. In each spatial bin, the 
continuum at the line was estimated using the flux in the 
adjacent energy band and the previously defined ratio. 

To complement this image approach, we present az- 
imuthally averaged radial profiles of the deconvolved line 
images. For the MOS1 and MOS2 cameras, the FWHM 
of the point spread function is 4.4" at 1500 eV and is 
slightly larger at higher energy (Aschenbach et al. 2000, 
Gondoin et al. 2000). 

We have used a standard Van Cittert deconvolution 
algorithm (Starck et al. 1998) to reconstruct the images 
within typically 20 iterations. 

3.1. Large scale structure 

The 15 A Fexvn line image is shown in Fig. 2, overlaid 
with the Si XIII He a contours. Overall, the spatial dis- 
tribution is very similar, although Fig. 2 shows that the 
image, as well as the radial profile, in the Sixin line is 
probably slightly peaked exterior to the Fexvn line im- 
age. The good correlation between the emission from iron 
and silicon implies that the outer layers of the supernova 
have been at least partially mixed. This is predicted by 
numerical hydrodynamical modeling of the global X-ray 
spectrum of Tycho (Itoh et al. 1988; Brinkmann et al. 
1989). 

In Fig. 3, we show again the 15 A Fexvn line image, 
overlaid with the contours in the continuum subtracted 
Fe Ka line image. There is an overall agreement between 
the line images, with spherical emission brightest in the 
northwest, faintest in the east and intermediate in the 
south sectors as for Sixin. Both Fe lines exhibit strong 
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Fig. 1. Top: EPIC spectra integrated over the whole rem- 
nant obtained with the MOS 1 (green), MOS 2 (red) and 
PN (black) camera. Bottom: RGS 1 spectrum of a limited 
region of Tycho (see Fig. 2) with order 1 and order 2. The 
spatial extent of the SNR dominates the line profiles. 



emission in the southern part of the southeast knots, and 
a faint contribution to the northern part for Fexvn. The 
comparison of the radial profiles and of the images demon- 
strates that the spatial distribution of the Fe Ka emission 
is broader, and peaks at smaller radius, than the Fexvn 
emission, as was first suggested by Hwang and Gotthclf 
(1997), based on ASCA deconvolved narrow band images. 

The fact that the Sixin He a and Fexvn images 
correlate well overall, whilst the Fe L and Fe K lines 
of iron are clearly displaced, favors the interpretation 
of a spatial structure in the temperature and ionization 
timescale in the shocked ejecta. This is indeed what is ex- 
pected from hydrodynamical evolution models of young 
supernova remnants (see for example Chevalier 1982 and 
Decourchelle & Ballet 1994). Although a spatial segrega- 
tion of the elements is possible, depending on the degree 
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Fig. 2. Top: Contours of the Sixm triplet blend (1670- 
2000 eV) superimposed on the EPIC MOS image in the 
Fexvn (A = 15 A, 775-855 eV) line. The RGS field of view 
is overlaid and the dispersion axis is parallel to the solid 
white lines. The Right Ascension is in decimal degrees. 
Bottom: Comparison of the azimuthally averaged radial 
profiles of the deconvolved images in the Si K (solid line) 
and Fexvn lines (dashed line). 
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Fig. 3. Top: Contours of the Fe K line (6200-6600 eV) im- 
age superimposed on the EPIC MOS image in the Fexvn 
line. Fields CI, C2 and C3 are shown with red solid lines. 
Bottom: Comparison of the radial profiles of the decon- 
volved images in the Fe K (solid line) and Fexvn lines 
(dashed line). 



of mixing, a structure in temperature and/or ionization 
stage is in most cases unavoidable and depends on the 
initial density profiles of the supernova. 

The Fe K line corresponds to higher temperature ma- 
terial (at least twice higher) and lower ionization timescale 
than the Fe L and silicon emission (Hwang et al. 1998). 
This indicates that the temperature in the inner shocked 
ejecta is higher with a lower ionization age than the outer 
shocked ejecta region. Such profiles of temperature and 
ionization age are those qualitatively expected from the 
evolution of a supernova remnant of 1.4 Mq in a uniform 
interstellar medium, when the reverse shock is propagat- 
ing, from an outer power-law density profile, into the in- 
ner plateau region of the ejecta. For an explosion energy 
of 10 51 ergs and an ambient density of 0.35 amu cm~ 3 , 
the reverse shock, at the age of Tycho, has already en- 
tered the central plateau (Decourchelle 1994). The dense 
shocked ejecta region (colder and older), close to the con- 
tact discontinuity, is the result of the interaction of a 
power-law density ejecta profile with the ISM (Chevalier 



1982) at the beginning of the evolution. When the reverse 
shock reaches the central plateau, the density begins to de- 
crease in the ejecta, the temperature increases and the X- 
ray emissivity diminishes (see for example Fig. 6 in Band 
& Liang 1988). This is qualitatively what we observe in 
Tycho with the inner Fe K line in a region with fainter 
emission. 

Standard models of the Type la supernova explosion, 
such as the W7 model (Nomoto et al. 1984), give outer 
density profiles in the ejecta closer to an exponential de- 
cline than to a power law. However, the characteristic tem- 
perature profiles for exponential ejecta expanding into a 
constant density ambient medium are flatter than the tem- 
perature gradient required in Tycho, unless a circumstellar 
medium is invoked (Dwarkadas & Chevalier 1998). Note 
however that for stratified ejecta, the heavier element lay- 
ers would have a higher temperature because of the lower 
number of electrons for the same mass density. 
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Fig. 4. Top: MOS 1 and MOS2 Spectra of the three 
southeastern knots CI (red), C2 (black) and C3 (green). 
Bottom: Spectra of a strong continuum region (Knot A in 
red) and a line dominated region (Knot B in blue). The 
regions A and B are shown in Fig. 5. 



Table 1. Derived abundances (relative to the solar pho- 
tospheric values) in the knots 
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tween the knots, clear variations of the abundances of sil- 
icon, sulphur and iron, which are given in Table 1. 

The spectral analysis demonstrates that the different 
emission structures seen in the southeastern knots are pri- 
marily due to spatial variations of the relative abundances 
of silicon, sulphur and iron, rather than temperature and 
ionization timescale effects. The relative abundance of sul- 
phur to silicon is roughly constant in these three knots, 
while that of iron is varying. Knot CI has the largest sili- 
con and sulphur abundance, whereas knot C3 has almost 
no iron. In this particular knot C3, almost a factor of 10 
less iron has protruded into the silicon layer compared 
with the adjacent knot C2. However, the southeastern 
knots are the only place where such a relative variation 
of silicon to iron is unambigously observed. As shown in 
Fig. 2 and discussed in section 3, there is otherwise a good 
agreement between the images in Fe XVII and in Si xni for 
the whole remnant. These knots resemble to some extent 
the silicon rich shrapnels seen in Vela (Aschenbach et al. 
1995; Tsunemi et al. 1999). 



3.2. Small scale inhomogeneities 

A more careful inspection of the comparison between the 
silicon and iron images (see Fig. 2) shows that of the south- 
eastern knots (Vancura et al. 1995) the northern one is 
brightest in the Sixni line, whereas the southern one is 
brightest in the Fexvn line. 

Fig. 4 (top panel) shows the MOS 1 and MOS 2 spec- 
tra of the three southeastern knots, marked CI, C2 and 
C3 on Fig. 3. We have extracted spectra from these iden- 
tically sized regions and fit them with a two component 
non-equilibrium model (with variable abundances for the 
first component). A Gaussian line at 3.1 keV was added 
to the model to account for the argon emission line, which 
is missing in the non-equilibrium model VNEI within 
XSPEC (Arnaud 1996). 

The first component is at a relatively low temperature 
(Ti) and accounts for the line emission which is the focus 
of this discussion, whilst the second component assumes 
solar abundances, is at a high temperature and accounts 
for the high energy continuum. We assume no hydrogen 
and helium in the shocked SN la ejecta, and solar values 
for carbon and nitrogen. 

The three different spectra can be simultaneously fit- 
ted with the same temperature (kTi = 0.62 ± 0.02 keV), 
ionization timescale (1.2 x 10 11 s cm -3 ) and hydrogen col- 
umn density (0.83 x 10 22 cm -2 ) and a relatively good fit 
is achieved (\ 2 = 1183/721 dof). However, there are, be- 



4. High-energy continuum and radio emission 

The high-energy continuum image in the 4500-5800 eV 
band is shown in Fig. 5, overlaid with the Sixni emis- 
sion contours. The image in the continuum is more regular 
and extended than in the silicon line, which shows irreg- 
ular outer contours, which can be attributed to Rayleigh- 
Taylor instabilities at the contact discontinuity between 
the ejecta material and the ambient medium (see Fig. 6). 
The extent of the continuum emission corresponds well 
with the radio VLA image at 22 cm (Dickel et al. 1991) 
and gives the position of the outer blast wave. This is 
shown in Fig. 5, where the radio image is overlaid with the 
contours of the continuum emission. The regular spheri- 
cal emission of the continuum, unlike the emission line 
images, may indicate that the ambient medium surround- 
ing Tycho was fairly homogeneous, although it has been 
recently interacting on its eastern edge with Hi clouds 
(Reynoso et al. 1999). The more asymmetric shape in the 
line emission may be related to azimuthal variations of the 
heavy element distribution in the ejecta. 

The radial profiles of the X-ray continuum and radio 
emission are shown in the bottom panel of Fig. 5. The peak 
of the continuum emission is the outermost and is inter- 
preted as emission from the shocked interstellar medium. 
The radio emission extends as far as the continuum but the 
profile peaks at about the same radius as the Si emission 
(Fig. 2). For comparison the profile of the Fe K emission 
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Fig. 5. Top: High energy continuum image (4500-5800 eV) 
overlaid with Sixin contours. Fields A and B are shown 
with black solid lines. Middle: VLA 22 cm radio map 
overlaid with high energy continuum contours. Bottom: 
Comparison of the three corresponding radial profiles: ra- 
dio (dotted), high energy continuum (dashed) and Fe K 
line (solid line). 



is also shown and clearly peaks at a much smaller radius 
than the other components. 

In addition to a regular smooth profile, the continuum 
map shows striking bright knots on the eastern edge, as 
well as on the diametrically opposite side of the remnant. 
These features are lying inner to the continuum edge, they 



Fig. 6. Deconvolved silicon image overlaid with outer con- 
tour of the radio image from Dickel et al. (1991). 



have no clear counterpart in the line images, nor in the ra- 
dio. However, the brightest continuum knot (eastern edge) 
lies very close to the two brightest radio peaks, which ex- 
hibit the flattest spectral index of the remnant (Filament 
I, in Katz-Stone et al. 2000), while the southwestern con- 
tinuum knots lie along a radio structure, which also has 
a flat spectral index (Filament A, in Katz-Stone et al. 
2000). The eastern knot has the lowest radio expansion 
in Tycho (Reynoso et al. 1997, see Hughes 2000 for the 
X-ray expansion rate), and is close to one of the brightest 
optical filaments (Kamper & Van den Bergh 1978; Smith 
et al. 1991), associated with infrared emission at 10.7-12 
micron (Douvion et al. 2000). A recent Hi 21 cm study of 
the environs of Tycho has found a strong H I absorption 
feature, corresponding to a cloud of density 160-325 cm -3 
on the eastern edge (Reynoso et al. 1999). 

The bright continuum knot lies at the periphery of the 
line emission image, whilst being well inside the continuum 
emission. The comparison of the deconvolved silicon image 
with the radio image shows that the eastern side is the only 
place where the radio and silicon boundaries overlap and 
where such a strong X-ray continuum is seen (see cross 
mark in Fig. 6). The two brightest radio knots coincide 
with similar features in the silicon image. This indicates 
that they are related to ejecta material, with structures 
characteristic of Rayleigh- Taylor instabilities as argued by 
Velasquez et al. (1998), based on VLA radio images with 
a resolution of 1". 

Jun & Norman (1996) have shown that Rayleigh- 
Taylor instabilities amplify strongly the ambient local 
magnetic field around the dense features of the instabil- 
ity by stretching, winding and compressing the magnetic 
field lines. Assuming an external magnetic field tangential 
to the shock front, their simulations account for the ob- 
served radial vector polarization observed in the remnant 
(Dickel et al. 1991; Reynoso et al. 1997) and predict an am- 
plification of the magnetic field in front of such Rayleigh- 
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Taylor features, which would explain the fact that these 
structures are seen both in radio (amplified synchrotron 
emission around these instabilities) and in X-rays (ther- 
mal emission of the dense knots). 

The strong continuum knot might be the result of the 
interaction of the eastern edge with H I clouds, leading 
either to higher densities and magnetic field strengths 
(hence stronger nonthermal emission), or to a reflected 
shock in front of the cloud, which heats and compresses 
a second time the ISM (hence leading to an increase 
in the thermal bremsstrahlung emission at high energy). 
However, no H I features have been observed in the dia- 
metrically opposite edge, where a strong X-ray continuum 
is also observed. 

The bipolar symmetry of the regions emitting this 
strong X-ray continuum may be an indication of particle 
acceleration with an external magnetic field tangential to 
the shock front in these regions as in the case of SN 1006 
(Reynolds 1996, 1998). 

Fig. 4 (bottom panel) shows the MOS 1 and MOS 2 
spectra of Knot A and Knot B (highlighted in Fig. 5). The 
spectrum of Knot A features a strong high energy contin- 
uum, whereas the spectrum of Knot B is typical of the 
brightest silicon dominated regions. We have used a single 
non-equilibrium model to fit the data (VNEI, XSPEC). 
The silicon and sulphur abundances in Knot A are close 
to solar, but all other elements have much lower values, 
in particular iron (see Table 1). The temperature is more 
than a factor 2 higher in Knot A (kT = 2.3 ± 0.3 keV) for 
an ionization timescale (n e t ~ 1.4 ± 0.4 x 10 10 s cm' 3 ) 
three times shorter than in Knot B (kT = 0.9 ±0.1 keV, 
n c t ~ 4.7 ± 0.3 x 10 10 s cm" 3 ). 

Both the narrow band images and the spectral analy- 
sis indicate that the continuum dominated emission most 
probably arises from the shocked ambient medium. We 
cannot, however, unambiguously identify the nature of 
this emission. The signature of a thermal plasma in this 
spectrum comes primarily from the silicon and sulphur 
lines, although they can also correspond to spatial con- 
tamination (due to the PSF) from the much brighter sur- 
rounding ejecta emission. A deeper systematic analysis of 
all the strong continuum knots is required to determine 
precisely their origin, which is beyond the scope of this 
paper. 

5. Conclusion 

We present the first results obtained from the XMM- 
Newton observation of the Tycho supernova remnant. 

- The good correlation between the images in the Sixin 
K line and Fe xvn L line implies that some fraction of the 
inner iron layer has been well mixed with the outer silicon 
layer. 

- The fact that the Fe K line emission peaks distinctly at 
a smaller radius than the Fe L line emission is a clear indi- 
cation that there is a spatial structure in the temperature 
in the ejecta, with higher temperature towards the reverse 
shock. This is coherent with the propagation of the reverse 



shock in the inner plateau of the ejecta. 

- The silicon image correlates well at small scales with the 
radio image, and probably marks the contact discontinu- 
ity, distorted by Rayleigh- Taylor instabilities. 

- The high-energy continuum, which is regular overall, 
peaks just behind the shock front, and is attributed to 
the emission from the shocked ambient medium. 

- In addition to this regular continuum emission, we have 
discovered particularly bright and hard continuum knots 
on the eastern and western edges. 
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